Radially Excited States of IP Charmonia and X(3872) 
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The excited states of charmonia are numerically investigated in quenched lattice QCD with im- 
proved gauge and Wilson fermion actions formulated on anisotropic lattices. Through a constrained 
curve fitting algorithm, the masses of the first excited states in ++ , 1 ++ , and l" 1 channels are de- 
termined to be 3.825(88), 3.853(57), and 3.858(70) GeV, respectively. Furthormore, a node structure 
is also observed in the Bethe-Salpeter amplitude of the 1 ++ first excited state. These observations 
indicate that X(3872) could be the first radial excitation of Xrf- 

PACS numbers: 12.38.Gc, 14.40.Lb, 11.15.Ha 
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The narrow charmonium-like X(3872) (with width F < 
2.3 MeV) was first observed by Belle in the exclusive 
decay B ± -> K±X -> K^+Tt-J/tp Q, and has been 
confirmed by CDFII 0], D0 and BaBar Q in three 
decay and two production channels. Even though its J PC 
quantum numbers have not been finally established, the 
present experimental data strongly favor that it is a 1 ++ 
state [5] for the following reasons. First, the decay mode 
V(3872) — > jJ/ip observed recently by Belle requires the 
charge conjugation of Y(3872) to be positive Q. The 
possibility with J PC = ++ and h can be ruled out 
based on the angular correlations in the 7r+7r~ J/ip sys- 
tem [l| and that of 2 ^ or 1 ^ is also strongly disfa- 
vored according to the dipion mass distribution [6| . Sec- 
ond, the bound r(e+e-)Br(X tt+tt-J/^) < 10 eV 
at 90% C.L., obtained with the data collected by BES 
at y/s = 4.03 GeV Q, implies that it is unlikely a 1~~ 
vector state. 

Since the discovery of A" (3872), there have been many 
theoretical studies. Given its quantum number J PC = 
1 ++ , a natural assignment of the state is the first radial 
excitation of IP charmonium state Xd- However, there 
are two main difficulties for this interpretation. One is 
its tiny decay width relative to other charmonium states, 
another is that it lies roughly 100 MeV lower in mass 
than the prediction of the quark model [|| . These diffi- 
culties motivate many non-charmonium explanations of 
V(3872), such as hybrids glueballs [lOfl . diquark clus- 
ters [llj]j and molecular states 

Although the non-relativistic quark model is successful 
for heavy quark bound states, it is known that relativistic 
effects can be important for charmonia , since the charm 
quark is not heavy enough. In particular, these effects, 
which are not taken into account in the model, can be 
more important for higher excited states. It is therefore 
more desirable to study charmonia with a relativistic lat- 
tice QCD formalism, which includes all the relativistic 
effects. In contrast to the study of the ground states of 
charmonia, such as the IS* and IP states, their excited 



states have not been investigated as much in the formal- 
ism of lattice QCD, even though they are more interest- 
ing in the present era when many new heavy mesons of 
open-charm and closed-charm are observed. The major 
obstacle is that the extraction of the excited states re- 
mains a challenge in Monte Carlo simulations. In this 
work, we investigate the relevant charmonium spectra in 
quenched lattice QCD and focus on the derivation of the 
first (radially) excited states through the sequential em- 
pirical Bayes method (SEB) [l3[ advocated by xQCD col- 
laboration, which is in the spirit of the constrained curve 
fitting algorithm and has been successfully applied to the 
study of nucleon excited states [l4| and pentaquarks [l5| . 



TABLE I: The input parameters for the calculation. Values 
for the coupling f3, anisotropy £, the lattice spacing a s , lattice 
size, and the number of measurements are listed. 
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We use the quenched approximation in this study. The 
gauge configurations are generated by the tadpole im- 
proved gauge action [l6[ on anisotropic lattices with the 
temporal lattice much finer than the spatial lattice, say, 
£ = a s /at 1, where a s and at are the spatial and 
temporal lattice spacing, respectively. Each configura- 
tion is separated by 2000 heat-bath updating sweeps to 
avoid the autocorrelation. The much finer lattice in the 
temporal direction gives a high resolution to hadron cor- 
relation functions, such that masses of heavy particles 
can be tackled on relatively coarse lattices. The relevant 
input parameters are listed in Table [J , where a s 's are 
determined from r^ 1 = 410(20) MeV. 

For fermions we use the tadpole improved clover ac- 
tion for anisotropic lattices The parameters in the 
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action are tuned carefully by requiring that the physical 
dispersion relations of vector and pseudoscalar mesons 
are correctly reproduced at each bare quark mass [Hj]. 
The conventional interpolation operators ijjTip are used 
for meson states with different gamma matrices T for 
the specific spin-parity quantum numbers. After the 
Coulomb gauge fixing, the wall-source quark propagators 
are calculated under the anti-periodical boundary condi- 
tion and are used to construct the correlation functions 
(the point-source correlation functions are found to be 
very noisy in ++ , 1 ++ , and 1 H channels). The bare 
charm quark masses at different (3 are determined by the 
physical mass of J/ip mj/^ = 3.097 GeV. 
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FIG. 1: The SEB fitting procedure of the vector charmo- 
nium correlation function at /3 = 2.8. The upper panel shows 
the fitted masses using fitting models with 1-4 mass terms. 
The lower panel shows the X* /d.o.f with the change of fitting 
range. 

A typical correlation function takes the function form, 

C(t)=J2Wie- Mit , (I) 
»=i 

where Wt and Mi are the spectral weight and the mass 
of the i-th state, respectively. If one fits the correla- 
tion function using this function form with multiple mass 
terms directly through the conventional maximal likeli- 
hood method, usually one gets a poor result due to the 
complicated parameter space and the result depends cru- 
cially on the the choice of the initial values of the param- 
eters. A constrained curve fitting algorithm can help if 
one can obtain reasonable priors of the parameters that 
will be fitted (l9j . In this work, we apply the sequential 
empirical Bayes method (SEB) in the data analysis. We 
will outline the main ideas of this method here and more 
details can be found in Ref. [13j. 



It is obvious in Eq. (jTJ) that the contribution of the 
lowest state dominates the correlation function when 
t is large enough, because the relative contribution of 
higher states to the ground state damps exponentially as 
cx exp(-AMjit), where AMji = Mj - Mj. Intuitively, 
if W\ and M\ of the first state are correctly derived in 
a time region [ti,t max ], one can treat them as priors for 
the ground state and then do the two-mass-term fit in a 
larger time range [t? , t max ] with ti < f i . The criterion 
for including the second mass term is the observation 
of a sharp jump of x 2 /n.o./, which is a signal that a 
single exponential cannot describe the data well, when 
t\ is decreased further. This procedure is repeated by 
adding more states until the data points are exhausted. 
This is the basic fitting procedure of SEB method. Gen- 
erally speaking, the last state can not be taken as the 
realistic one because it includes almost all the contami- 
nations from higher states. Figure [T] illustrates the SEB 
fitting procedure in the vector channel. The upper panel 
shows the fitted masses using 1-4 mass terms in the fit 
function, while the lower panel shows the x 2 per degree 
of freedom (% 2 /d.o.f). We perform one-mass fit in the 
time range \t\,t max \ = [50,68], and then add the second 
mass term at t — 49 where a sharp increase of x 2 /d.o.f 
is observed. In the two-mass fit, the fitted mass of the 
second state becomes more and more stable when de- 
creasing t2, and the x 2 /d.o.f does not change much all 
the way down to t = 30 = t<z- At t — 29 x 1 /d.o.f begins 
to climb up and we add the third state from there to do 
a three-mass fit. Similarly, we add the fourth state from 
t = 14 down to t — 4 and find that the four-mass fit 
can model the correlation function very well in the whole 
time range [t m in,t max ] = [4,68] with a x 2 /d.o.f < 1. 
One can see from the figure that not only the mass of 
the ground state keeps constant, but also the masses of 
the second and even the third state are also very stable 
after the fourth state is included in the fitting model. The 
best fit results of the masses for the three lowest states 
are 0.4367(2), 0.517(3), and 0.577(10), which correspond 
to the physical masses 3099(2)MeV, 3669(22)MeV, and 
4084(90)MeV, respectively, given a s = 0.139fm. These 
results, which are in good agreement with the experi- 
mental values, illustrate the efficacy of the SEB method. 

Now we come to the analysis of the correlation func- 
tions in ++ , f ++ , and l" 1 channels. It is found that 
these correlation functions are much noisier than those 
in vector (V) and pseudoscalar (PS) channels, such that 
we can only carry out maximally three-mass-term fits to 
them in the practical data analysis. Figure shows the 
fitting procedure for the 1 ++ charmonium at (3 = 2.8. 
The upper panel is similar to that of Figure [TJ but shows 
the three-mass-term fit, where one can see that the sec- 
ond state is very stable after the inclusion of the third 
state in the fitting model. The lower panel illustrates the 
effective mass plateau m e g = \n C (t) / In C(t + 1), where 
the data points are the simulation results while the curve 
is the three-mass fitting function with the best-fit param- 
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FIG. 2: The fitting procedure in 1 ++ channel. The upper 
panel is similar to that of Figure [1] but shows the three-mass- 
term fit. The lower panel illustrates the effective mass plateau 
m e fl = In C(i)/ In C(t + f), where the data points are the 
simulations results while the curves are the fitting function 
plotted using the best fitted parameters. 



eters. The good agreement of the curve and the data 
points manifests the quality of the three-mass fit. The 
fitting procedure for the ++ and l" 1 channels is similar. 

We also carry out the same calculations at (3 = 2.4, 
2.6 with smaller statistics (200 configurations each) to 
consider the continuum extrapolation. Due to the larger 
lattice spacings at these two /3's as shown in Table [H 
the largest time for the fit windows in parity-positive 
channels is roughly t max = 20 and 30 respectively. For 
both j3, the correlation functions can be well fitted by 
only two mass terms in the full time windows and the 
third mass terms are found to be marginal. This implies 
that the fitted first excited states may have significant 
contaminations from higher states. The continuum limits 
are obtained by the linear extrapolation in terms of a 2 s . 
Table |TT] lists the best-fit masses of the ground and the 
first excited states at the three /3's and their continuum 
extrapolations. The extrapolation values of the -P-wave 
ground states are in agreement with previous works [20l . 
[21| but lower than the experiment values. The results of 
the extrapolated masses of the first excited states should 
be taken with caution because the masses at the two 
small /3's can have large systematic errors as mentioned 
above. 

The masses of the excited masses obtained at j3 = 2.8 
are relatively more reliable and can be treated as approxi- 
mations of their continuum value with the awareness that 
the uncertainties from the finite lattice spacing are not 
properly tackled. Comparing with the predictions by the 
non-relativistic quark model [8| (the last colume of Ta- 



TABLE II: Best-fit masses of the ground and the first excited 
states of parity-positive charmonia at different j3. All errors 
are statistical. The errors of the continuum limit values are 
from the linear extrapolation in terms of a 2 s . 



2.4 



2.6 



2.8 



cont. 



++ (lP) 3.466(19) 3.437(21) 3.410(18) 3.376(6) 
1 ++ (1P) 3.520(19) 3.507(17) 3.477(14) 3.453(6) 
1+-(1P) 3.515(17) 3.505(17) 3.488(13) 3.472(8) 



++ (2P) 3.816(57) 
1++(2P) 3.937(56) 
1+~(2P) 3.955(53) 



3.865(90) 
3.887(54) 
3.964(71) 



3.825(88) 
3.853(57) 
3.858(70) 



3.857(55) 
3.800(2) 
3.835(88) 



BGSI8] 



3.424 
3.505 
3.516 



3.852 
3.925 
3.934 



blc|TT|, the mass of the 1 ++ first excited state at (3 = 2.8 is 
lower than quark model expectation but consistent with 
the experimental value of X(3872). Even though our re- 
sult is not decisive to the assignment that A" (3872) be the 
first radial excitation of Xci due to the large statistical 
error, it indicates that the possibility that X(3872) is a 
conventional charmonium state can not be simply ruled 
out. 

In order to see if the extracted masses are really those 
of the first excited states, we also investigate the Bcthc- 
Salpeter amplitudes of charmonium states at (3 — 2.8. In 
the Coulomb gauge, we split the sink operator into two 
parts, with each quark field residing on different spatial 
sites, namely, Or{x,y) — iJj(x)Tip(y), where T represents 
the various gamma matrices corresponding to specific 
J PC quantum numbers and x — (x, t) and y = (y,t). 
In the practical calculation, all the two-point functions 
with |x— y| < 6\/3a s are calculated. It is found that the 
two-point functions exhibit a spherical symmetry with 
respect to the relative displacement between the quark 
and anti-quark fields, r = x — y. Therefore the two- 
point functions with the same spatial separation |r| are 
averaged to increase the statistics. 

The two-point functions calculated with the 
Coulomb wall sources are actually Cr(r, t) = 
E(0|O r (x,t;x + r,i))O^(0)|0), where O r ,w{0) = 



£^(y,o)iY(z,o). 



After the integration over x, we get 



C r (r,t) = ]r-^-<0|O r (0,0;r,0)|n) x 



2M 

n 

£>|O{,(y,0;z,0)|0>e 

£*»(»•) 



-M„t 



-M„t 



(n=l,2,...), (2) 



where M n is the mass of the n-th state and (r) the BS 
amplitude (up to an irrelevant prefactor) we would like 
to extract. 

The r-dependent BS amplitudes $ n (r) have direct con- 
nections with the non-relativistic wave functions of heavy 
quarkonia [13, HH . Specifically, for the P-wave charmo- 
nia, $„(r) is related approximately to the radial wave 
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FIG. 3: $ n (r)'s for the ground (n — 1) and the excited state 
(n = 2) in the 1 ++ channel. The mini-panel shows the en- 
larged plot around r/a — 8. 

function R n (r) as 

$„(r) ~ Ci£(r), (3) 

where R' n {r) is the r-derivative of R n (r) and C a fac- 
tor irrelevant to the discussion here. In this work, <E>i(r) 
and $2( r ) are obtained from a constrained-curve-fitting 
algorithm by assuming the fitted masses of the ground 
and the first excited states at r = as known parame- 
ters, say, priors. Plotted in Figure [3] are $ n (r)'s of the 
ground (n — 1) and the first excited state (n — 2) in the 
1 ++ channel. It is obviously seen in the figure that <E>i(r) 
has one node at r/a ~ 4.8 and $2( r ) has two nodes at 
r/a ~ 2.5 and r/a ~ 8. These behaviors are qualitatively 
in agreement with the theoretical anticipation described 



by Eqn.Q if both states correspond to the IP and 2P 
states of 1 ++ charmonium: the unique node of $i(r) 
corresponds to the unique maximum of IP radial wave 
function, the two nodes of $2(f) manifest the two ex- 
trema, and necessarily a radial node in-between, of the 
2P radial wave function. 

To summarize, we have carried out a quenched lattice 
study of the excited states of charmonia on anisotropic 
lattices. In our study, a relativistic formalism of lattice 
QCD is employed and hence the relativistic effects are 
included. Using SEB, which is in the spirit of the con- 
strained curve fitting method, we can derive reliably the 
masses of the ground states and the first excited states 
of charmonia in the ++ , 1 ++ , and l" 1 channels. We 
obtain a mass of 3853(57) MeV for the first excited state 
of the 1 ++ charmonium, which is lower than the conven- 
tional quark model predictions but consistent with the 
measured mass of A" (3872) within statistical errors (sys- 
tematic errors from the quenched approximation have not 
been considered yet). We have also observed a radial 
node of the wave function of the 1 ++ first excited state. 
All these results support the possibility that X(3872) be 
the first radially excited state of \ci ■ 
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